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We show how a candidate mean-field amplitude can be constructed from the exact wave function 
of an externally forced A^-Boson system. The construction makes use of subsidiary [N — l)-particle 
states which are propagated in time in addition to the true A-particle state, but does not involve 
spontaneous breaking of the 17(1) symmetry associated with particle number conservation. Pro¬ 
vided the flow in Fock space possesses a property which we call maximum stiffness, or t-coherence, 
the candidate amplitude actually satisfies the time-dependent Gross-Pitaevskii equation, and then 
serves as macroscopic wave function of the forced A-particle system. The general procedure is il¬ 
lustrated in detail by numerical calculations performed for the model of a driven bosonic Josephson 
junction, which allows one to keep track of all contributions which usually are subject to uncon¬ 
trolled assumptions. These calculations indicate that macroscopic wave functions can persist even 
under conditions of strong forcing, but are rapidly destroyed upon entering a regime of chaotic dy¬ 
namics. Our results provide a foundation for future attempts to manipulate, and actively control, 
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I. INTRODUCTION AND OVERVIEW 


This paper concerns the time evolution of a system of 
N interacting identical spinless Bose particles of mass to, 
where N is assumed to be large. At each moment t in 
time it is described by a Schrodinger wave function 'F 
which is symmetric in all its spatial arguments, 


^'(ri,...,rAr;t) = T'(rp(i),..., rp(jv); t) (1) 


for each permutation P of the indices 1,... ,N. Its dy¬ 
namics are governed by an V-particle Schrodinger equa¬ 
tion 




( 2 ) 


N 




i=i 




T'(ri,...,rAr;t) 


where 

hi{r;t) =-^A + V{r;t) (3) 

denotes a single-particle Hamiltonian containing an ex¬ 
ternal potential V which may depend explicitly on time, 
so as to exert some controlling influence on the system, 
and the interparticle interaction U naturally is symmet¬ 
ric. 


U{r,r') = U{r',r) . (4) 

We will assume that U can be replaced by an effective 
contact interaction with strength g, 


( 5 ) 


which means that we restrict ourselves to the low-energy 
regime [T]. 

The central question addressed in this work is under 
what conditions an approximate description of the dy¬ 
namics in terms of a mean-held amplitude $ evolving 
according to the time-dependent Gross-Pitaevskii equa¬ 
tion [2H9] is viable, and, if so, how that amplitude is 
related to the actual wave function df. Since $ = <i)(r;t) 
depends on only one spatial coordinate r, it contains 
much less information than the many-body wave func¬ 
tion 0 could possibly carry. Therefore, a reduction of 
the full iV-particle dynamics to the mean-field level with¬ 
out essential loss of information is feasible only if already 
the 7V-particle wave function itself is relatively simple, 
or ordered, so that the mean-field amplitude, if it exists, 
may be regarded as an order parameter of the system m- 
According to common knowledge deriving from London’s 
theory of superfluids m the mean-field amplitude should 
correspond to a “macroscopic wave function”, that is, to 
a macroscopically occupied single-particle orbital. This is 
fully in line with the preservation of information: In the 
ideal case where all N particles occupy the same orbital, 
no information is lost if $(r;t) equals that orbital. But 
this leads to further questions when the system is sub¬ 
jected to time-dependent forcing: Even if the initial state 
is given exactly by an iV-fold occupied single-particle or¬ 
bital, this initial order might be destroyed by the ex¬ 
ternal force, and the time-dependent many-particle state 
may become arbitrarily complicated, no longer admit¬ 
ting a mean-field description. On the other hand, one 
may formally solve the Gross-Pitaevskii equation with 
an arbitrary type of external forcing incorporated into 
the single-particle Hamiltonian hi. Hence, there must be 
some sort of indicator which tells one whether or not the 
solution to the time-dependent Gross-Pitaevskii equation 


Ues{r,r') = gS{r- r') 
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can actually serve as a macroscopic wave function of the 
forced A^-particle system. Here we explore this intuitive 
idea in mathematical terms. 

As a guideline for the following deliberations, and to 
state the subject as clearly as possible, we briefly sketch a 
popular pseudo-derivation of the Gross-Pitaevskii equa¬ 
tion. To this end, we introduce the usual Fock space an¬ 
nihilation and creation operators ij{r) and ijj^r) which 
obey the canonical Bose commutation relations 


tp{r),'ip{r') = 0 

= 0 

il){r),ip\r') = 5(r — r') 


( 6 ) 


so that the Fock-space many-body Hamiltonian takes the 
form [12] 

H{t) = j (^r'ip\r)hi{r-,t)'ij}{r) (7) 

-I- - y d^r J d^r' ip\r)tj}'^ {r')U{r,r')'tl}{r')ip{r) . 

Here we use the “hat”-symbol to designate operators act¬ 
ing in Fock space, in contrast to single-particle operators 
such as the one given by Eq. (§. 

With the help of the unitary operator U{t, to) which ef¬ 
fectuates the time evolution from some initial moment to 
to t and thus obeys the equation 


ifi-^U{t,to) = H{t)U{t,to) , 


( 8 ) 


the field operator is transformed to the Heisenberg 
picture through the familiar prescription 

'^H{r;t) = U\t,to)${r)U{t,to) , (9) 

leading to an equation of motion of the form 

^{r\t),H(t) U{t,to) . (10) 

Working out the commutator appearing on the right- 
hand side, one obtains 




( 11 ) 


+ / '$l^{r'-t)U{r,r')il)^{r'-t)'ip^{r;t) ; 


inserting the contact potential (§ , one is left with 
ih^$ii{r-t) = hi{r-,t)i)ii{r-,t) 


( 12 ) 


Now comes a major assumption: The mean-field ampli¬ 
tude $ is supposed to be given by the expectation value 
of the Heisenberg field operator. 


VN^{r;t) := 


(13) 


This definition requires clarification. If it is taken liter¬ 
ally, each Wparticle state |4')Ar, no matter how “simple”, 
can only produce a left-hand side equal to zero, since the 
(A^—l)-particle state i^H{r;t)\d')N is orthogonal to |'I')jv 
in Fock space. Hence, when working with Eq. ( |13| ) one 
supposes that the superselection rule for the total particle 
number N is somehow violated, and the state with 
respect to which the expectation values (13) are taken 


is some coherent superposition of Wparticle states 
with different N, 


I'k) — ^ aAr|'I')Ar 


(14) 


N 


This corresponds to the idea that the U{1) symmetry 
associated with particle number conservation is sponta¬ 
neously broken. While that concept may be helpful for 
simplifying formal calculations, it remains an auxiliary 
device not rooted on the grounds of the actually given 
A^-particle system. This deficiency has led to several 
profound discussions in the literature, and to the devel¬ 
opment of number-conserving approaches |SH3 EMU- 
The unreflected use of the notion of spontaneous sym¬ 
metry breaking has been severely challenged by Leggett, 
who criticizes that there are no circumstances in which 


Eq. {' 14 ) is the physically correct description of the sys¬ 


tem, and maintains that the definition (131 is “liable to 


generate pseudoproblems, and is best avoided” |7]. Thus, 
one of the goals of the present study is to provide an 
alternative, logically more transparent definition of the 
mean-field amplitude for an A^-Boson system. In contrast 
to customary analysis, our approach does not involve a 
splitting of the field operator into a condensate part and 
a remainder. Instead, we will first introduce a candidate 
mean-field amplitude by taking matrix elements of the 
field operator with the physical A^-particle state and cer¬ 
tain auxiliary (N — l)-particle states, and then provide 
a criterion which decides whether or not that candidate 
is physically meaningful, so that it qualifies as a macro¬ 
scopic wave function. 

Still staying within the framework of spontaneously 
broken symmetry, an appealing possibility to satisfy 
Eq. (13) seems to arise when Idt) is a coherent state, 
which, by definition, is an eigenstate of the annihilation 
operator: 


^iiir;to)\dJ) = VN(^{r;to)\d>) , 


(15) 


or, upon return to the Schrodinger picture, 

V((r)|4'(to)) = //V4>(r;to)|^'(to)) ; (16) 

in this respect such many-body coherent states are 
analogs of the familiar standard harmonic-oscillator 
coherent states already introduced in the early days 
of quantum mechanics [m in]. Intuitively speaking, 
Eq. (15) means that jT) is “simple” in the sense that 
it remains unchanged if one particle is removed. How¬ 
ever, if 14') is an eigenstate of ^h(J';1o) at one particu¬ 
lar moment to, it will not, in general, be one at others. 
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Therefore, we are confronted with the necessity to admit 
more general concepts of “coherence”. 

Following the traditional route for the time being 
and taking the symmetry-broken expectation value (131, 
Eq. @ immediately yields 


ih—'/N^{r;t) = hi{r;t)'/N^{r;t) 

+ 9 


(17) 


This leads to another pertinent problem; Of course one 
could write down the equation of motion for the third mo¬ 
ment appearing here, but this would involve still higher 
moments, and so on, leading to an infinite chain of equa¬ 
tions expressing an infinite momentum hierarchy |18Lll9j . 
The Gross-Pitaevskii equation results from the assump¬ 
tion that this chain can be terminated at the lowest pos¬ 
sible level, so that the expectation value of the triple 
operator product is replaced by the product of the ex¬ 
pectation values of the individual operators: 


{tfjR{r;t) 


(18) 

Here the “~^”-sign is meant to indicate that this closure 
hypothesis (18) generally constitutes an uncontrolled 


approximation; note, however, that it would be exact 


for a coherent state (15). Accepting this factorization, 
Eq.([^ finally turns into the celebrated Gross-Pitaevskii 
equation [21M IT9U22] . namely. 


= hi{r;t)^{r;t)-\-Ng\^{r;t)\‘^^{r;t) . (19) 
at 

Since one has the identity 

/'dV(^lj(r;t)^H(T’;t)) = Af (20) 


mean-field amplitude is normalized to unity, 

y'dV|$(r;t)|2 = l, 


for any t, one deduces from the definition (|13|) that the 

( 21 ) 


now assuming factorization of the expectation value (20). 


Besides providing an alternative definition of the mean- 
field amplitude not based on spontaneous breach of sym¬ 
metry, questioning and monitoring the accuracy of the 
required closure hypothesis constitutes a main objective 
of this work. 

We proceed as follows: In Sec. [TT] we outline a ten¬ 
tative construction process of a candidate mean-field 
amplitude associated with an Wparticle state |'I'(t))Ar. 
This construction rests on the introduction of subsidiary 
{N — l)-particle states |'I'(r|t)) 7 v-i which are obtained 
from |T(t))jv by the annihilation of one particle; the dy¬ 
namics of these states are governed by the same Hamil¬ 
tonian as those of the physical state |'I'(t))w- Our ap¬ 
proach is similar in spirit to the introduction of a conden¬ 
sate wave function by Lifshitz and Pitaevskii in Ref. [53] , 


but places particular emphasis on the way the initially 
“neighboring” states |d>(t))Ar and |'I'(r|t))jv_i evolve in 
time: Do they, in a suitable sense, remain close to each 
other? The subsequent derivation of the candidate’s 
equation of motion then makes clear under what con¬ 


ditions the Gross-Pitaevskii equation (19) is satisfied, so 


that the construction yields a bona fide macroscopic wave 
function. If, on the other hand, these conditions are not 
met, the candidate bears no physical meaning and has 
to be discarded. These considerations are not meant as 
a substitute for rigorous mathematical analysis [T3l422j . 
but rather as an elaboration of the salient features which 
make the Gross-Pitaevskii equation work. In particular, 
we isolate a notion of coherence which does not corre¬ 
spond, in the first place, to a single-particle orbital being 
macroscopically occupied, but which places a certain de¬ 
mand on the flow in Fock space instead, requiring its 
stiffness. 


The general scheme is applied in Sec. |IH| to a nu¬ 
merically tractable model system, describing a driven 
Josephson junction. The model is presented in Sec. [Ill A] 
and the corresponding equations of motion are derived 
in Sec. [Ill B[ now paying particular attention to keep¬ 
ing all error terms. The results of selected numerical 
calculations which support our reasoning in considerable 
detail are discussed at length in Sec. |HIG[ Besides be¬ 


ing paradigmatically simple, the model has the further 
pleasant feature that it lends itself to a semiclassical in¬ 
terpretation, allowing one to relate the A^-particle dy¬ 
namics to those of a driven nonlinear pendulum. As laid 
out in Sec. EnDl this semiclassical view leads to impor¬ 
tant qualitative insights which may guide the analysis 
of more complicated, experimentally accessible systems 
in the future. Thus, the model calculations collected in 
this Sec. jnT] complement recent studies of periodically 
(5-kicked Bose-Einstein condensates which aim at explor¬ 
ing nonlinear dynamics in a quantum many-body con¬ 
text [24H29] : here we emphasize the potential usefulness 
of a sinusoidal force with a slowly varying envelope. 

We finally phrase our main findings, and draw some 
general conclusions, in Sec. EYl 


II. CONSTRUCTIVE APPROACH TO THE 
MEAN-FIELD AMPLITUDE 


We start by constructing the V-particle Fock-space 
state vector |'I'(t))Ar corresponding to some given wave 
function Building on a vacuum state |vac) character¬ 
ized by 


and 


(vac I vac) = I 


'i/’(r)|vac) = 0 for each r 


( 22 ) 


(23) 
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we have 

I^W)jv = -^Jd^ri ... Jd^rN'i'iri,... ,rN;t) 

X ijj\ri) .. .tp\rN)\veLc) . (24) 


another moment ti, whereas the Gross-Pitaevskii equa¬ 


Repeatedly using the identity 

N 




i=i 


N 


= ^(5(r'-rfc) (25) 


k=l 




which follows directly from Eqs. and ( |^ , one con¬ 
firms that this state is properly normalized, 

^(vE-(t)|'I-(t))^ = l, (26) 

provided the Schrodinger wave function is, 

yd^ri... J |4'(ri,..., t)|2 = 1 . (27) 


Schrodinger equation 


This state vector (|24|) evolves in time according to the 

(28) 


= H{t)\d>{t))N 

with the many-body Hamiltonian ([7). 

The system’s one-particle reduced density matrix, as 
given by 

Q{r,r']t) = jv(4'(t)|^'*'(r)^(r')|4'(t))7v , (29) 

then takes the explicit form 

g{r,r';t) = N Jd^r 2 ... Jd^VN (30) 

X 4'*(r,r2,...,rAr;t)4'(r',r2,...,rAr;f) . 

It might appear natural to introduce a time-dependent 
macroscopic wave function on the basis of the Penrose- 
Onsager criterion for Bose-Einstein condensation in an 
interacting iV-particle system [3D]: At each moment to 
one can, in principle, perform a spectral decomposition 

g{r,r';to) ='^n,{to)<f*{r,to)<f^{r',to) (31) 


of this density matrix (29) into a complete orthonor¬ 


mal set of single-particle eigenfunctions (pi{r,to)- If then 
only one of the eigenvalues, say no(to), is on the order 
of N, while all other ones are at most of order 0(1), 
one has a simple Bose-Einstein condensate at to, and the 
eigenfunction ipo{r,to) can be regarded as the instanta¬ 
neous macroscopic wave function. However, according 
to this construction the phase of (po(r,to) remains un¬ 
determined, and unrelated to the phase of (po{r,ti) at 


tion (19) for $(r;t) constitutes an initial-value problem 


which does not leave such a freedom. Hence, even in the 
most ideal case in which no (to) = N for all to, it is not 
possible to reconstruct the full time-dependent macro¬ 
scopic wave function from the one-particle reduced den¬ 
sity matrix, even if one knows the exact many-body state 
at all times, because the phase of the macroscopic wave 
function evidently requires additional consideration. 

Therefore, here we take a different route. In situations 
where a mean-field amplitude $ can be meaningfully de¬ 
fined, the most obvious requirement to be placed on such 
an object is that its absolute square, multiplied by N, 
be equal to the diagonal elements of the density matrix, 
that is, to the expectation values of the density operator, 
at least to a good approximation: 

7V|4>(r;t)p = Ar(4'(t)|'^'l'(r)^(r)|4'(t))Ar . (32) 

Given an initial state |4'(to))Af at some arbitrary initial 
moment to, it is easy to satisfy this requirement on a for¬ 
mal level at that moment to by invoking the normalized 
(N — l)-particle states 

|^(r|to))^-i := , (33) 

where 7 (r) is a phase that can be suitably adjusted later 
on, and we have 

||V^(r)|4/(to))^f 

= N J d^r2 ■ ■ ■ J d^rN \'i>(r,r2, ■ ■. ,rN;to)\‘^ (34) 


according to Eq. (30), as corresponding to the particle 


density at the position r. Thus, it may happen that the 
denominator appearing in the definition (33) vanishes at 


certain points r; such points either are exempted from 
the analysis, or require special treatment. 

It is quite essential not to be misled here: |'I'(r|to)) n-i 
must not be regarded as some sort of wave function de¬ 
pending on r; rather, for each r it is a state vector in the 
(N — l)-particle sector of the Fock space. We will mark 
such subsidiary (N — l)-particle states by the “tilde”- 
symbol in the following. 

Introducing the projection operator 

Pto = \^ir\to))N-i Ar-i(5(T-|io)l , (35) 

one obviously has 

AT (4'(<o) IV''*'(r)V^(r) |4'(to)) JV 
= N{dt(to)\^^Hr)Ptoi’ir)\d>{to))N (36) 


by the very definition ( |33[ ); as a consequence, we now 
may set 

VN^{r,to) ■■= 7v-i($(r'|to)|^(T’)|^(to))Af • (37) 
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Evidently the projector (35) does not depend on the 


phase 7 (r) of the subsidiary states, but this phase al¬ 
lows one to tune the phase of ^{r,to). 

This tentative construction (|3^ still does not mean 
that the expression $(r,to) actually is a macroscopic 
wave function at the moment to- The above prescrip¬ 
tion could be applied to any initial state |'I'(<o))Ar, re¬ 
gardless of whether or not this state admits a mean-field 
description. At this point, a suitable counterpart of the 
Penrose-Onsager eigenvalue criterion which distinguishes 
a condensate is still lacking. 

The usefulness of this construction becomes evident 
only when it is extended to later times: Besides the evolv¬ 
ing physical A^-particle state 


|'I'(f))Ar = U{t,to) |'I'(to))7V 


(38) 


we also have a family, parametrized by r, of auxiliary 
{N — l)-particle states which are obtained by subjecting 
the subsidiary initial states (33) to the evolution gener¬ 


ated by the very same Hamiltonian, giving 

|$(r|t))Ar-i = U{t,to) |$(r|to))jv-i • (39) 

We now define a candidate mean-field amplitude $ for 
all t > to by taking, in literal extension of the initial-time 
definition (37), the matrix elements of the field operator 
with these evolving states (p^ and (39): 


\/]V$(r;t) := Ar_i(^'(r|t)|^(r)|4'(t))Ar . (40) 

Note that, in contrast to the “symmetry-breaking” def¬ 
inition (13), this still tentative definition is free of arbi¬ 


trariness. The right-hand side of Eq. (40) can in principle 
be evaluated once the initial state |'k(t ^ n is given; there 
is no freedom comparable to the choice of the coefficients 
qn in Eq. (14). However, it is not yet clear at this stage 


whether the quantity d* defined in this manner actually 
conforms to the central requirement (32) for t > to, and 


whether or not it satisfies the Gross-Pitaevskii equation. 
Obviously, the answers to these questions depend to a 
large extent on |'I'(tQ)) 7 v, but also on the particular time- 
dependence imposed on the Hamiltonian. On the other 
hand, it is easy to see that our construction does fulfill its 
purpose in the ideal case where all N particles initially 
occupy the same orbital ip: Then the initial many-body 
wave function is given by a Hartree product state. 


N 

5'<p(T’i,...,rjv;to) = ; 

1=1 


(41) 


in this case one expects ^{r,to) to be equal to (p{r,to). 
Indeed, following the general prescription (24) the 


state |4>,^(to))Ar associated in Fock space with this prod¬ 
uct wave function (41) acquires the form 




1 




d^rip{r;to)'fHr 


N 


I vac) (42) 


and thus obeys the quasi-eigenvalue equation 
^{r)\'if^{to))N = VN(p{r;to)\'i>^{to))N-i , 


(43) 


implying that the one-particle reduced density ma¬ 
trix (29) here becomes 


p(r, r'-to) = Nif* (r; to)(p{r'- to) ■ 


(44) 


This reflects precisely the Penrose-Onsager criterion for 
Bose-Einstein condensation [30]: In the case of a pure 
condensate the one-particle reduced density matrix (31) 


equals the N-fold of the projection operator onto the N- 
fold occupied single-particle orbital. Moreover, Eq. (43) 
immediately yields 


i){r)\d>^{to))N ‘f{r;to) 


||V'(r)|d'^(to))7v|| Mr;to)\ 


|'I'<^(to))Ar-l 


(45) 


Hence, if one chooses the phase 7 (r) as the negative phase 
of the single-particle orbital p{r, to), so that, for instance, 
7 (r) = 0 when ip{r,to) is real and positive, Eq. (33) here 
reduces to 




3(’'l^o))Ar-l — l^i^(^o))iV-l 


(46) 


In conjunction with the quasi-eigenvalue equation (43) 
and our definition ( |40| ) of the mean-field amplitude, this 
furnishes the expected initial identity 


d>(r,to) = <p(r,to) ■ 


(47) 


What may be even more important, Eq. (46) also dras¬ 
tically simplifies the entire construction process for all 


subsequent times t > to, because the family (39) reduces 


to a single member not depending on r, and for all times 
the dependence of $ on r in Eq. (40) is inflicted solely 
by the field operator ^{r). 

The Hartree product wave function (41) embodies 
the idea of “simplicity” in a particularly straightforward 
manner. If the many-body wave function carries merely 
the information content of a single-particle orbital, no in¬ 
formation is lost when employing that orbital as a mean- 
field amplitude, in the spirit of a macroscopic wave func¬ 
tion. For the sake of easy nomenclature we will refer to 
the associated Fock-space states (42) as Wparticle co¬ 


herent — in short: N-coherent — in the following. This 
wording appears to be somewhat at odds with the quasi¬ 
eigenvalue equation (43), which is reminiscent of the ac¬ 


tion of an annihilation operator on an harmonic-oscillator 
eigenstate, or on a typical Fock state, in marked distinc¬ 
tion to the proper eigenvalue equation (16) characterizing 


a conventional coherent state. But it is this type of co¬ 
herence which is implied by macroscopic occupation with 
a sharp number of particles. Needless to say, an initially 
iV-coherent state generally will not remain iV-coherent 
in the course of time. Thus, the question what property 
an iV-particle state must have in order to admit a Gross- 
Pitaevskii description of its time evolution still needs to 
be settled. 
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Returning to the general case not necessarily starting 
from an TV-coherent state (42), the equation of motion of 
our candidate mean-held amplitude (401 is given by 


iS—y]V$(r;t) = Ar-i(^'(r|t)| ^{r),H{t) |^'(t))jv . 

(48) 

In formal agreement with the previous Heisenberg equa¬ 
tion the commutator yields 


ip{r),H{t) = hi{r\t)^{r) (49) 

-I- [d^r' ijj\r')U{r,r')ijj{r')'ijj{r) , 


having exploited the symmetry Now the hrst term 
on the right-hand side leads to the matrix elements 

hi(r;t) Ar_i($(r|t)|'^(r)|4'(t))jv 
= hi{r-,t)'/N^{r\t) , (50) 


alternate between the (TV — 1)- and the TV-particle sector: 

J'ip\r')U{r,r')'tp(r')'ij){r) 

= Jd^r' U{r,r')^{r)^'^{r')ip{r') 

— jd?r'U{r,r')5{'i — r')'ijj{r') . (51) 

The hrst contribution appearing here is dealt with by 
replacing the actual interaction potential by the effective 
contact interaction ([^, giving the triple operator product 
{r)ijj{r). Taking the matrix elements required by 
the dehnition ( [4^ , we then write 

g N-i{^ir\t)\ipir)ip\r)ip{r)\d)(t))N 
5 Ar_i($(r|t)|^(r)|^'(T))jv 

X Ar(^'(T)|'^'l'(r)|$(r|T))Ar-i 

X Ar_i(T(r|T)|'0(r)|4'(t))Ar 

= 5TV3/2|<I>(r;T)|2$(r;t) . (52) 


where, in analogy to Eq. ( |18[ ), the -symbol indicates 
that this particular step can be correct for some spe¬ 
cial states, but will lead to uncontrolled errors for most 
others. That is, steps marked by this symbol determine 
about the quality of our construction: The magnitude 
of the error committed there decides whether or not it 
will result in an acceptable mean-held amplitude. In 


the present case (501 the single-particle Hamiltonian hi , 


when willfully drawn before the matrix element in order 
to recover the candidate mean held amplitude (40), also 


acts on the r-dependence of the auxiliary family (39) 


— unless there is no such dependence. At this point, 
the observation made below Eq. (47), deriving from the 


property (46) of TV-coherent states, leads to a noteworthy 
insight: 


If the initially given N-particle state is ex¬ 
actly N-coherent, in the sense specified by 
Eqs. (fl) and (4-2), it suffices to employ 
the corresponding (TV — l)-particle state (46) 


as sole subsidiary state, implying that the 
chain \5(f( is an equality for all later times 
t > to. 


Thus, one may reasonably hope that the error admit¬ 
ted when working with close-to-TV-coherent initial states 
remains tolerably small, but substantial further work 
seems required to turn this vague hope into a quanti¬ 
tative statement. 

The second term on the right-hand side of Eq. (49) 
confronts us with a different problem: Acting twice with 
the field operator on an TV-particle state would bring us 
to the (TV — 2)-particle sector of Eock space, over which 
we have no command. Therefore, with the help of the 
commutation relations we reorder the operators such 
that their successive application produces images which 


Note that appears once again: Here the projec¬ 

tion operators projecting onto the states |4'(r|T))and 
respectively, have been inserted between the 
field operators “as suitable”, but so far without deeper 
justihcation, thus providing a formal substitute for the 
standard closure assumption ( [I^ . We will investigate 
the validity of this step exemplarily in Sec. |HI[ 

The second contribution arising from Eq. ( |5I[ ) necessi¬ 
tates to evaluate the interaction potential U{r,r') at zero 
distance, r = r', and thus forbids us to employ the effec¬ 
tive substitute (§. This appears to be not more than a 
spurious complication, so that here we restrict ourselves 
to interaction potentials which remain finite at zero dis¬ 
tance, such as a simple repulsive step potential [3T] . Then 
U(r,r) ^ eg will be about proportional to the coupling 
constant g, with a factor c carrying the dimension of an 
inverse volume. Thus, after taking the required matrix 
elements we are left with 


eg jv-i(^'(r|t)|^(r)|4'(t))7v = cgy/N^{r-,t) . (53) 

Collecting the elements (50), ([^, and ([^, the basic 


equation of motion (48) adopts the form 


= hi(r;T)$(r;T) 


Ng (\^{r;t)f --^^^{r-t) . (54) 


Taking the limits TV —>■ oo and g —>■ 0 such that the 
product Ng remains constant, this actually becomes the 
Gross-Pitaevskii equation (19), confirming that our def¬ 


inition ( |40[ ) indeed has the potential to work as desired, 
still with the proviso that both the initial state |'I'(to))Ar 
and the Hamiltonian H(t) be such that the replacements 
indicated by in Eqs. (50) and (521 constitute at least 
good approximations. 
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However, to rederive the Gross-Pitaevskii equation was 
not what we wanted to achieve in the first place. Rather, 
we had set out to construct a mean-field amplitude 
such that its absolute square yields the exact expectation 
values of the density operator according to Eq. (321. Is 
this the case now? 

Evidently, in order to factorize the right-hand side of 
Eq. (32) exactly we should be writing 


AT (r)'0(r) l^'(f)) AT 

with projection operators 

Qt = |S(r|f))Ar-i Af-i(S(r|f)| 
constructed from the normalized states 

^(r)|^'(f))Ar 


(55) 


(56) 


i(r|f))jv-i = 


IIV'(r')|'k(i))iv|| 

ll^(r-)|'k(f))^|| 


(57) 


which are obtained by Hrst propagating |4'(fo ))n in time, 
and then annihilating a particle at r. Instead, Eq. (32) 


combined with the definition (40) implies that the success 


of our construction would be guaranteed if one had the 
identity 


AT (4'(f) (r)v)(r) |4'(f)) JV 

where, in continuation of Eq. ( |35[ ) , the operator 
Pt = |$(r|f))Ar-i Ar-i(5(r’|t)| 
projects onto the ray generated by 

||V^(r)|4^(fo))^|| 


(58) 


(59) 


(60) 


In contrast to the previous state (57), this is the state 


found by first removing a particle from |'I'(fo))Ar at 
and then propagating from fg to t. These two projection 
operators Pt and Qt are identical, meaning that Eq. (58) 


actually is an exact equality, if |S(r|f))Ar-i differs from 
|4>(r|f))Ar-i merely by a phase factor, and thus is pro¬ 
portional to |4'(r|f))Ar-i. This latter demand requires 


^(r)|4'(f))Ar = t) |4'(r|f))Ar-i 


(61) 


with some function r/(r; f); taking the scalar product with 
|$(r|f))Ar-i immediately yields 


(62) 


The resulting consistency condition 
'0(r){7(f,fo)|4'(to))Ar 

= , (63) 


or, in more compact form, 

V)(r)|4'(f))jv = V)/V$(r;f)|$(r|f))Ar-i , (64) 

constitutes an interesting generalization of the quasi¬ 
eigenvalue equation (43). Whereas that former rela¬ 


tion characterizes fV-particle coherence at one fixed mo¬ 


ment to only, its descendant (64) intrinsically incorpo¬ 


rates the system’s time evolution: Annihilating a parti¬ 
cle from the time-evolved state should produce the state 
obtained by annihilating prior to evolution; states having 
this property will be called t-coherent from here onwards. 

Thus, we may phrase the preceding considerations as 
follows: 


If the N-particle state under consideration is 
t-coherent, in the sense specified by Eqs. (63) 


or (64), the absolute square of the mean-field 
amplitude defined by Eq. equals the par¬ 
ticles’ density, obeying Eq. !(3S^. 


Obviously, the desired identity (64) selects those cases 


in which a good macroscopic wave function does exist. 
Even if it is satisfied initially to a good approximation, it 
may cease to be valid under the action of external forc¬ 
ing. The value of the above conceptions now rests in the 
recognition that they furnish a tool for monitoring the 
achievable performance of the Gross-Pitaevskii equation: 
What really matters is the degree of t-coherence, as ex¬ 
pressed by the extent to which the states (57) required 
for exact factorization equal the states (|60[) underlying 
our construction; this extent is quantified by the abso¬ 
lute value of the “direction cosine” EH 


R{r;t) = e‘'^(’'V-i(^'(T’|i)|S(r|<))Ar-i (65) 

^ ^i-y(r) N-i{^{r\t)\i((r)\'i/{t))N 

^ N{'If(.to)\f>^r)U'<{t,to)f({r)U{t,to)\-^{to))M 

||^(r)|vI/(to))^|| ||^(r)|4/(f))^|| 

Namely, if |i?(r;t)| = 1 the solution $(r;f) to the Gross- 
Pitaevskii equation does provide an accurate description 
of the system; if 0 < |i?(r;f)| < 1, it is bound to miss 
some aspects of the full iV-particle dynamics. In that 
sense, R{r;t) quantifies the “degree of mean-field ap- 
proximability” of the evolving system at time t, and at 
the position r. Dynamically speaking, the magnitude 
|i?(r;f)| evaluates the evolution of two trajectories in 
Pock space initially differing by one particle; if it does 
not matter whether one annihilates a particle first and 
then lets the system evolve over some interval of time, 
or evolves first and annihilates then, the flow in Fock 


7y(r;t) = •\/]V4)(r;f) . 
























space may be considered stiff. Therefore, we refer to the 
magnitude |i?(r;t)| as stiffness^ with maximum stiffness 
\R{r;t)\ = 1 indicating perfect t-coherence. Note that 
such maximum stiffness still does not mean that annihi¬ 
lation and evolution commute in the operator-theoretic 
sense, but rather that the system manages to hide the 
non-commutativity of these operations for the state con¬ 
sidered behind a phase factor which drops out of the pro¬ 
jectors involved. 


It is also of interest to observe that if Eq. (64) is satis¬ 
fied the second line of Eq. (65) yields 


R(r-1) = 


( 66 ) 


so that the phase of the solution to the Gross-Pitaevskii 
equation, which is relatively easy to obtain numerically, 
contains information about how |S(r|t))jv_i differs from 
|$(r|<))Ar_i, while both of these states designate the 
same ray in Fock space; this information quantifies the 
difference between the evolution of an N- and that of an 
{N— l)-particle state. On the one hand, this observation 
ties up with the well-known fact that in thermodynamic 
equilibrium, when the Hamiltonian is time-independent, 
the phase of the macroscopic wave function is given by 
the chemical potential, i.e., by the energy required to add 
one more particle to the system. On the other hand, it 
also clarifies why it is not possible to determine the phase 
of the macroscopic wave function from the one-particle 


reduced density matrix (31), even if the Penrose-Onsager 


criterion is optimally satisfied: That matrix draws its in¬ 
formation solely from the fV-particle sector of Fock space. 

It goes without saying that the theoretical construc¬ 
tion process outlined in the present section is not meant 
as a practical help for solving iV-particle evolution equa¬ 
tions; rather, it provides a conceptual framework within 
which time-dependent many-Boson systems can be ana¬ 
lyzed. The following study of a model which allows one 
to carry out all of the above steps at least numerically 
will confirm that the general terms coined here do in¬ 
deed capture essential features of iV-Boson dynamics in 
a nontrivial manner. 


III. CASE STUDY: THE DRIVEN TWO-MODE 
SYSTEM 


A. The model 

We now investigate in detail a system of N identi¬ 
cal spinless repulsively interacting Bose particles which 
are confined to two sites, labeled 1 and 2. Each pair of 
Bosons occupying the same site is assumed to invariably 
contribute the repulsion energy 2 hK to the total energy of 
the system, while intersite interaction is neglected. More¬ 
over, there is a tunneling contact between the two sites, 
with hopping matrix element Ml 12, so that B denotes the 
single-particle tunneling frequency. In terms of the oper¬ 
ators aj and which create and annihilate, respectively, 


a particle at the jth site, obeying the Bose commutation 
relations 




= 0 


a] >4 


= 0 




= 5 


'jk 


(67) 


with j,k = 1,2, the system then is described by the 
Josephson Hamiltonian [7] 


Hn=-- 


nn 


(^a\a2 + +hK (^a\a\aj^a^ + a\a\a2ar 


( 68 ) 

This may be regarded as a particular variant of a model 
introduced by Lipkin, Meshkov, and Click in order to test 
the validity of many-body approximation methods [53- 
155] - in a spirit quite similar to that of the present paper. 
It is realized approximately by a Bose-Einstein conden¬ 
sate loaded into an optical double-well potential, thus 
providing a “bosonic Josephson junction” |36j : in this 
context the model has met with substantial renewed in¬ 
terest pTZl - HH] . For ease of notation we write its Fock- 
space operators without “hat”-symbol. 


Obviously the general field operator ^{r) with its con¬ 
tinuous “index” r, as considered in the previous section, 
here corresponds to the pair oi, 02 . Likewise, the macro¬ 
scopic wave function <I>(r;t) now reduces to a set of two 
discrete mean-field amplitudes ci(t), C 2 (t), the absolute 
squares of which, multiplied by N, should give the num¬ 
bers of particles residing at the two sites. Introducing 
the dimensionless time variable r = Qt, their equations 
of motion, easily derived by the scheme sketched in the 
Introduction, take the form 


i^Ci(T) = -ic2(r)-F2a|ci(r)|2ci(r) 

i^C2(T) = -ici(r)-F 2a|c2(r)pC2(r) , (69) 

depending only on the scaled interaction strength 

Nk 


a = 


B 


(70) 


This system (691 constitutes the Gross-Pitaevskii equa¬ 


tion for the Lipkin-Meshkov-Glick model (68); it also de¬ 
scribes self-trapping on a polaronic dimer 113 iZ] 


It 


is well known that this system can be integrated ana¬ 
lytically in terms of Jacobian elliptic functions gaig. 
But such integrability clearly signals that the model (68) 
still is too simple for our present purposes: For test¬ 
ing the viability of a Gross-Pitaevskii-type approxima¬ 
tion to full V-particle evolution under typical conditions 
we require a system which exhibits more generic, non- 
integrable mean-field dynamics. For this reason we ex¬ 
tend the model by adding an explicitly time-dependent 
driving term: Respecting the general form of a “canoni¬ 
cal Josephson Hamiltonian” [7], we study the driven two¬ 
mode system 




( 71 ) 
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where 


Hiit) = sin(a;t) {a\ai — a^a.^ (72) 


models a bias applied with frequency uj between the 
two sites, while its amplitude h^(t) may vary in time. 
Even with fixed amplitude /i(t) = this driven bosonic 
Josephson junction possesses chaotic mean-field solu¬ 
tions 


Adapting the construction (24), a general A^-particle 
state of the model can be written in the form 

■ ■ ■ “L l^ac) . 

Jiv —1,2 

(73) 

Because of bosonic symmetry the coefficients '0ii, 
coincide for all permutations of their indices. Writing 
= '>Pn{t) when n of the N indices 
equal 1 (so that the other N — n indices equal 2), we then 
have 

\m). = ;7^En^.w(4)”(4)""”ivac) 

^ ■ n—O ^ 

= \n, N - n) , (74) 

n=0 ^ 


where 


N-n 


|n, N — n) = 


(<_( 4 )_ 

a/uI — n)\ 


I vac) 


(75) 


denotes the Fock state with n particles on site 1 and, 
consequently, N — n particles on site 2. Thus, a general 
Wparticle state is specified by -I- 1 time-dependent 
complex coefficients ipnit). 

An A^-particle coherent state which corresponds, in the 
sense of Eq. (42), to an A^-fold occupied single-particle 


orbital is characterized by the requirement that each cre¬ 
ation operator be accompanied by the same amplitude. 
Denoting the single-particle basis of the two-site system 
as {(, 91 , (^ 2 }, and parametrizing a general single-particle 
orbital tp in this basis by means of two angles 9, (j) in the 
form 


ip = cos{9/2)ipi -I- sin( 0 / 2 )e “^(/92 , (76) 

thus omitting an irrelevant overall phase, the coherent- 
state coefficients become 

= cos”(0/2) sin^-"(0/2)e‘(^-")‘^ ; (77) 

the A^-particle coherent states of the two-site model are 
therefore written as 

= EL) cos"(0/2) 

xsin^-"(6l/2)e*(^-”)‘^|n,Af-n) . (78) 


These states, which emerge here as special cases of the 
A^-coherent states (42), play a significant role in quan- 

“5'C/(2)- 
. In the 


turn optics, in which field they are referred to as 
coherent states” or “atomic coherent states” [S 3 ] 
present context their most important property lies in the 
fact that they satisfy the quasi-eigenvalue equations 


ai| 6 ',(/))Ar = '/Ncos{e/2)\9,(j))N-i 

a^\e,(l))N = VNsm{e/2)e^^\9,^)N-i , (79) 


in formal correspondence with the general Eq. (43). 


It is easy to see that the best A^-coherent approxima¬ 
tion to the ground state of the undriven Josephson junc¬ 
tion ( 68 ) is provided by the symmetric state 


k/ 2 , 0 )iv = 


1 

2 W 2 


N 

E 

n—O 


1/2 


,N — n) 


(80) 


giving 


n{7t/2,0\Ho\7t/2,0)n 

Nhn 


la/ 1 

2 2 N 


(81) 


with the scaled interaction strength (70). However 


this A^-coherent state ( |80| ) is an eigenstate of the _M2- 
proportional “hopping” 


term of the Hamiltonian ( | 68 | ), 
but not of its AK-proportional interaction term, which 


would be diagonalized by the Fock states (75) instead. 


Therefore, the A^-coherent state (80) is a good approxi¬ 
mation to the true ground state only as long as a <C 1 ; 
for stronger interaction the true ground state must be 
more sharply peaked around n = N/2. Because one has 


A^y/^ 

nj 


jvk) 


(82) 


for |n —A^/2| «C N/2, a good ansatz for the exact ground 
state is 


/o \ 1/4 N 

^exp(^-"f{n-N/2)‘^'^\n,N-n) 

(83) 

with variational parameter 7 , supposing N ^ 1 . An 
elementary calculation then yields 

N{'l\Ho\-f)N ^ 1^1 7 1 A 

Nhn 2 \ N 2 2 -fN'^ J 

a / 2 1 \ 

^2 \ N^ 7A12 J 

+ 0 {N-^) , (84) 


where the neglected terms of order 0{N do not de¬ 
pend on 7 . Hence, one deduces the optimal parameter 


7 opt — 


V2a + 1 
N 


(85) 
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0.45 0.50 0.55 

n/N 


FIG. 1: (Color online) Comparison of the exact ground state 
of the two-mode model (681 for a = 0.95 and N = 1000 (dots) 
with the corresponding variational state (8^ with optimal 
width parameter (851 (full line), and with the A^-coherent 
state (801 (dashed). Plotted are the squares of the respective 
expansion coefficients ip„ in the Fock-state basis. 


which properly reduces to 7opt « 1/fV when a <C 1, as 
expected from the approximation (82), but grows with 


increasing interaction strength so as to suppress fluctu¬ 
ations. Figure compares this variational ground state 
for a = 0.95 and N = 1000 to the numerically computed 
exact one, confirming the quality of the variational ap¬ 
proach. More generally, one has 


I 2(2a +1)1/4 

Ar(7r/2, 0|7opt)Af — -— 

V 2a -I- 1 -f 1 


( 86 ) 


for large N, implying that the ground state of the Joseph- 
son junction (68) is quite different from the TV-coherent 


state (80) for large a, and does not approach it in the 
limit iV —>■ oo when this limit is taken such that a re¬ 
mains constant. This will become important for assessing 
the numerical studies reported in Sec. |III C[ 


B. Equations of motion 

The goal now is to explicitly carry through the con¬ 
struction process outlined in Sec. |ITj to monitor its ac¬ 
curacy, and to keep track of the errors committed, for 
the driven two-mode system We are given the 

Schrodinger equation 

ifi^|vk(t))=7V(t)|vI/(t)) (87) 

for Wparticle Fock-space states |5'(t)) which are writ¬ 
ten without particle number index here, and we wish to 
construct mean-field amplitudes Cj{t) such that 

N\c,{t)\^ = (88) 


for j = 1,2, in compliance with the central require¬ 
ment (32). Starting from an initial TV-particle state 
I Ik (to)) prepared at time t = tQ, we follow the pre¬ 
scription (33) and construct the two subsidiary (N — 1)- 
particle states 


|vk,(to)) := 


|a,l'k(to))ll 


(89) 


where the phases can be chosen at will, and we are 
assuming that both sites are initially occupied, so that 
the denominators do not vanish. These states obviously 
permit the factorizations 

(fl'(to)laja^lfl'(to)) 

= («'(/o)|a]|$,(to))(5,(/o)|a,|fl'(to)) (90) 

and therefore allow one to define the candidate mean-field 
amplitudes at the initial moment: 


VNcj(to) := (fi'j(to)|aJ«'(to)) ■ 


(91) 


Then not only the actual initial state |'I'(to)) is subjected 
to the time evolution generated by as expressed by 

the original Schrodinger equation (87), but also the two 
subsidiary states (89), giving rise to two further evolution 
equations 




(92) 


Here a decisive feature of the Fock-space formalism is 
exploited: Regardless of whether it acts on the A^-particle 
sector or on the (TV — l)-particle sector of Fock space, 
the Hamiltonian JI(t) is the same. Taking the solutions 
to these equations ([^ and ((Ml), the candidate mean- 


field amplitudes then are defined for all times t > to in 
accordance with Eq. (40): 

VTVc,(t):=($,(t)|a^.|4-(t)) . (93) 

Their equations of motion are obtained from the identity 

m^($,(t)|a,.|vl/(t)) = ($,(t)| [a^,H{t)] |vk(t)). (94) 
Evaluating the commutator 

= — —a2+‘2HKa\aiai+Hfj,{t)sm{ujt)ai (95) 
and taking the required matrix elements, one is led to 
ih^VNci{t) = -^VTVc2(t)/i(t) 

-|-2TiK('ki(t)|a|a]^aj^ l'k(t)) 

sm{ujt)VNci{t) , (96) 

where we have introduced the quantity 
((ki(t)|a2|4'(t)) , 


fi(t) = 


(4'2(t)|a2|4'(t)) 


(97) 
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an analogous, still exact equation holds for C 2 {t). Ob- 
serve the role of the ratio (97): When operating on 


Eq. (95) with l'Pi(t)) from the left, and with |'I'(<)) from 
the right, the tunneling term proportional to 02 yields 
its numerator, whereas the desired amplitude 02 ( 1 ) is 
given by its denominator. Wilfully setting /i(t) = 1, 
and doing the same with the ratio / 2 (t) appearing in 
the equation for C 2 (t), corresponds to the uncontrolled 
step marked by in the general Eq. (50); retain¬ 


ing both fi{t) and / 2 (t) in Eq. (96) and its C2(t) 


counterpart therefore renders these equations exact, and 
allows one to control one type of error accepted in the 
usual Gross-Pitaevskii approach. In the special case that 
the initial state |'k(to)) is given by an atomic coherent 
state (78), the quasi-eigenvalue equations (79) ensure 


that |5'i(to)) = l^ 2 (to)) = having set 71 = 0 , 

72 = —(/f, and hence one deduces |d>i(t)) = |'l> 2 (t)) for all 
times t > to- Therefore, in this special case one actually 
has fi{t) = / 2 (t) = 1 , so that the above error does not 
occur. However, in all other cases it needs to be consid¬ 
ered. 


To proceed with Eq. (96), we now adapt the general 


steps (51) and (52) for processing the triple operator 


products: 

{^j{t)\a]a Jail'S {t)) 

= 

- (^,(t)|a^.|vk(t))(vk(t)|a]|$,(t))($,(t)|a^.|vl/(t)) 

-VNcj{t) 

= . (98) 

The error introduced when assuming the factorization 
following the “^’’-symbol is given by the difference 

A,(t) = ($,(t)|a^.a]a^.|vk(t)) - . (99) 

Hence, we have the exact evolution equation 


id , , 1 , , „ , . 


-f 2 ; 


,Nk 




N 


Ai(t) 


-^sin(u;t)ci(t) , 


Ar3/2 

( 100 ) 


together with its counterpart for C 2 {t). Finally, taking 
the limit N —>■ 00 such that 


Nk 

~n 


= a = const. 


( 101 ) 


we obtain the Gross-Pitaevskii equation for the driven 
two-mode system 0 . with the proviso that both 


fjit) 


( 102 ) 


and 


A,-(t) 

Ar3/2 


0 


(103) 


in this limit, for j = 1,2: Again employing the dimen¬ 
sionless time variable t = we then have 

i(^ci(T) = -ic2(r)-f 2 a|ci(r)pci(r) 




sin(^r) ci(r). 


i:^C2(T) = -ici(r)-f 2Q!|c2(r)pC2(r) 


dr 


2 

M('r) „. 

n 


sin(^r) C 2 (t) . (104) 


The key q uest ion, o f cou rse, is under what conditions 
the limits ( 102 ) and (103) actually are adopted, and how 


relevant these limits are when N is still finite, and kept 
fixed. 

A further question deriving from the discussion in 
Sec. [h] is to what extent the solution to this system (104) 
does comply with the basic requirement ( 88 ). In analogy 


with Eqs. (55) and (56) we now have the exact identities 


(4/(t)|aJa,.14/(1)) 

= (4/(l)|at|S,(l))(S-(l)|a,|4/(l)) 


(105) 


with 


oW) = 






a^C/(l,lo)|^'(^o)) 

l|a,|4/(t))|| 


(106) 




whereas we require 

(4/(1 

= (4/(l)|a]|4/,(l))(4/,(l)|a^.|4/(l)) 
C/(l,lo)a^|4'(lo)) 


with 


|4/,(l))=e'^/- 


,l'P(to))ll 


(107) 


(108) 


a phase factor, c 


The actuall y ful filled Eq. ( |105| ) is compatible with the 
desired Eq. (107) if |Sj(l)) differs from |4/j(l)) merely by 


emanding a relation 

aj|4'(l)) = ; 


(109) 


the definition (93) of the candidate mean-field amplitude 


then immediately yields 

Vj{t) = . 


( 110 ) 
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Thus, we are led to a set of two consistency equations 
analogous to Eq. (641 which embody the property of per¬ 
fect t-coherence: 




( 111 ) 


this is the property which guarantees the validity of the 
basic Eq. (88) underlying the entire construction process. 


Eollowing Eq. (65), the degree to which t-coherence 


actually does prevail, that is, the stiffness of Fock-space 
flow in the vicinity of the initial state, is quantified by 
the scalar products 

R,[t) = e'^^($,(t)|S-(t)) 


= e‘^^-- 




lsl'k(to))|| 




.( 112 ) 


If both |Sj(t)) and |4'j(t)) represent the same ray in Fock 


space, as is implied by the consistency condition (111), 


then \Rj(t)\ = 1; monitoring the absolute value of this 
scalar product (112) therefore allows one to assess, on 


the basis of the exact fV-particle solutions |4'(t)) and 
the {N — l)-particle solutions |4'j(t)), the possible accu¬ 
racy of a Gross-Pitaevskii-treatment. M oreov er, if close- 
to-perfect t-coherence is given and Eq. (Ill) is satisfied 
more or less exactly, then one has 






(113) 


so that the phase of Rj (t) equals the phase of the respec¬ 
tive mean-field amplitude, up to . Read in the reverse 
direction, this means that if a Gross-Pitaevskii approach 
is viable, then the phase of the easily obtainable mean- 
field amplitudes contains valuable information about the 
exact many-particle dynamics, more precisely, on the dif¬ 
ference between the evolution of N and (iV — 1) particles. 


C. Numerical studies 



t/T 


FIG. 2: (Color online) Scaled population imbalance {Jz}/N 
(blue line), closure error |Ai|/N^^^ (red line), and sti ffness 
I Ri I (black line) for the undriven Josephson junction (681 with 
a = 0.5 and N = 1000. The initial state is the Fock state 
|N, 0), which is N-coherent. Here and in the following figures, 
time t is measured in multiples of T = 27r/fl. 



FIG. 3: (Color online) As Fig. but with N = 10000. 


Since the dynamics of the driven Josephson junc¬ 
tion (0- when occupied with N particles, takes place in 
a merely (iV-l-l)-dimensional complex Hilbert space, solv¬ 
ing the equations of motion for moderately large N does 
not pose severe computational difficulties. For ensur¬ 
ing high numerical accuracy we employ a variable-order 
Adams PECE algorithm, of the type originally elabo¬ 
rated by Shampine and Gordon [54] . routinely reaching 
particle numbers N on the order of a few thousand. 

The diagnostic tools now at our disposal are the 
ratios (97), the errors of closure (99), and the stiff¬ 
nesses (112). While the first and the second of these 


tools allow one to trace sources of deviations from Gross- 
Pitaevskii behavior, the third one evaluates the degree of 


t-coherence of the evolving many-body wave function, ad¬ 
mitting a truthful Gross-Pitaevskii description only when 
|i?,(t)| = l. 

As preliminary examples for the performance of these 


tools we monitor the dynamics of the junction (68) in 
the absence of the drive 


As initial condition we 
select |4'(0)) = |A1,0), so that all particles occupy site 1 
at to = 0. Then |4'i(0)) = |A^ — 1,0) is given by Eq. (89), 


with 7 i = 0. Moreover, since the initial Fock state | A^, 0) 
equals an Al-coherent state (78) with 0 = 0 and arbitrary 
(/}, we may set |'I' 2 (t)) = |'I'i(t)) for t > 0. This implies 
fi{t) = f 2 {t) = 1, SO that all deviations from mean-field 
dynamics are solely due to the closure errors. Figure 
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depicts the scaled population imbalance 

{Jz){t)/N = (4'(t)|alai - 4a2|«'(t))/(27V) 


0.5 


(114) 


for N = 1000 and scaled interaction strength a = 0.5, to¬ 
gether with the scaled error of closure | Ai(t)|/fV^/^ and 
the stiffness |i?i(t)|, with time t being measured in mul¬ 
tiples of T = 2tt/Q. The quantities | A 2 (t)|/fV^/^ and 
|i? 2 (t)| behave quite similar to their respective counter¬ 
parts. One observes the familiar collapse of the oscillat¬ 
ing population imbalance, caused by dephasing due to 
the finite particle number, which is to be followed by re¬ 
vivals at later times, when the components of the wave 
function rephase miMlESlES]. The most significant er¬ 
ror of closure occurs during the collapse stage, leading to 
a stiffness which decreases to zero in an oscillating man¬ 
ner. When increasing the particle number to iV = 10000, 
while keeping a = 0.5 fixed so that, in accordance with 
Eq. (101), the interparticle interaction strength Hk is re¬ 


fe; 


-0.5. 



10 


t/T 


20 


30 


FIG. 4: (Color online) Scaled population imbalance {Jz)/N 
for the undriven Josephson junction (681 with a = 0.5 for N = 


duced by a factor of 10, the collapse proceeds more slowly, 
as plotted in Fig. The comparison of the A-particle 
dynamics for both N = 1000 and N = 10000 with the 
prediction 


1000 (black line) and N = 10000 (blue line), in comparison 
with the prediction of the Gross-Pitaevskii equation (69) (red 


line). In each case, all particles are initially located at site 1. 


/ 7 \ — f\ I /o flic;'; ratios (97). Figure shows the dynamics for 

( 2 )mf( ) — (^|ci( )| |C 2 ( )\ j / ( ) Q, _ Q_g 5 comparatively small particle number N = 


of the Gross-Piaevskii equation ( [6^ displayed in Fig. 
provides convincing evidence for convergence in the limit 
iV —>■ oo, when a is held constant: Given any moment 
ti > 0, the particle number N can be increased such that 
the true populations of the two sites, calculated from the 
iV-particle Schrddinger equation, differ by an arbitrary 
small amount from the corresponding mean-field data in 
the entire interval from to = 0 to ti in the example con¬ 
sidered here. But as pointed out in Sec. Ill A| this sit¬ 
uation is exceptionally simple insofar as the mean-field 
equation of motion (69) is integrable. 


a = 0.95 and comparatively small particle number N = 
100 in response to a pulse with maximum scaled ampli¬ 
tude /imax/^2 = 0.9. Here the error of closure remains 
fairly small until a few cycles before the pulse’s enve¬ 
lope reaches its maximum. During that same interval one 
finds both /i(t) « 1 and / 2 (t) « 1 to good accuracy, and, 
most importantly, both stiffnesses |i?i(t)| and \R 2 {t)\ re- 


For assessing the conceptually far more difficult, but 
generic cases of non-integrable mean-field dynamics we 
now turn to the driven Josephson junction ( [7T| ), and 
choose a Gaussian envelope function 


= MmaxSXp ( 


(116) 


so that the drive (72) models a single pulse with maxi¬ 
mum amplitude S/Xmax, carrier frequency ui, and width a. 
For all following calculations we set w/D = 1.0 and 
a/T = 10.0, where again T = 27r/r2. Thus, the driving 
frequency is at resonance with the single-particle tun¬ 
neling frequency, and the driving amplitude increases 
smoothly from almost zero to its maximum within about 
30 driving cycles, then decreasing within another 30 cy¬ 
cles. 

As initial condition for large negative times we now 
select the ground state of the undriven junction (68), 
which, according to Sec. |III A] differs from an A-coherent 
state even for arbitrarily large N. Hence, now 14/2(0) 
is different from |'I'i(0), and we have to keep track 



FIG. 5: (Color online) Dynamics of the driven Josephson 
junction ( |71[ ) with a = 0.95 and N = 100 under the in¬ 
fluence of a Gaussian pulse ( |116[ ) with maximum amplitude 
Mmax/D = 0.9, scaled frequency u/Q. = 1.0, and width 
a/T = 10.0. Shown are the scaled population imbalance 
{Jz)/N (blue line), the closure error |Ai|/A^^^ (red line), 
the real part Re(/i) of the ratio (97) (orange line), and the 
stiffness |i?i| (black line). The initial state is the ground state 
of the undriven junction (1681). 
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FIG. 6: (Color online) Comparison of the population imbal¬ 
ance for N = 100 under the action of the pulse monitored 
in Fig. 1^ (blue line) to the corresponding prediction of the 
Gross-Pitaevskii equation (1041 (red line). 


FIG. 8: (Color online) Real and imaginary part of the ratio 
/i for the pulse monitored in Fig. 


main close to unity before they drop to values fluctuating 
around about 0.5 in the second half of the pulse. This di¬ 
agnosis forecasts that the Gross-Pitaevskii equation will 
describe the pulsed iV-particle dynamics quite faithfully 
until about the pulse’s maximum, and then become un¬ 
reliable. This prediction is fully borne out by Fig. 
which compares the true population imbalance (114) for 
this pulse to the mean-held imbalance (115), as computed 
from Eq. (104): As long as \Rj{t)\ « 1, the simple Gross- 
Pitaevskii equation does a remarkably good job. 


A far more impressive example for the possible accu¬ 
racy of the Gross-Pitaevskii equation is given in Fig. 
Here the particle number has been increased from N = 
100 to N = 1000, while all other parameters are the same 
as in Fig.[^ again implying a reduction of the interaction 



FIG. 7: (Color online) As Fig. bnt with N — 1000. 


strength Hk so as to keep a constant. Throughout this 
pulse the closure error remains negligible, and both ratios 
fj{t) stay close to unity, as witnessed on a hne scale by 
Fig.i This results in almost perfect t-coherence of the 
evolution, even though the initial state is not fully N- 
coherent, guaranteeing excellent mean-held approxima- 
bility of the entire pulse dynamics. Since the maximum 
scaled amplitude — 0.9 is by no means small, 

as is evident from the quite signihcant response of the 
population imbalance, this hnding is far from trivial. 

According to Eq. (113) the phase of the scalar prod- 



FIG. 9: (Color online) Comparison of the phase of the mean- 
field amplitude Ci(t) (red lii^ with that of Ri{t) (black) for 
the pulse monitored in Fig. [A for which |i7i(t)| = 1 to good 
accuracy. The insets magnify the evolution of the phases at 
the beginning of the pulse, and around its maximum, where in 
view of Fig. 1^ the mean-field approximation is most critical. 
At the end of the pulse, the phases differ by 0.210 tt. 
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FIG. 10: (Color online) As Fig. but with N = 1000 and FIG. 12: (Color online) As Fig. but with N = 10000 and 
/Tmax/H = 1.5. /imax/fl = 1.5, Confirming the “sudden death” of the mean 

field already observed in Fig. |10| 


uct Rj{t) should determine the phase of the mean-field 
amplitude Cj(t) in case of perfect t-coherence, when 
\Rj{t)\ = 1. This is confirmed in Fig. for the pulse 
studied in Fig. which meets the above requirement to 
good accuracy: With 71 = 0, the phase of ci(t) practi¬ 
cally coincides with that of Ri{t), growing linearly with 
time. 

An opposite paradigm is captured by Fig. again 
for N — 1000, but now the maximum scaled amplitude 
has been increased to ^max/f^ = 1-5. Here the response 
of the A^-particle system is characterized by a stiffness 
which remains close to perfect almost up to the pulse’s 
maximum, but then rapidly drops to values close to zero. 



as a consequence of a suddenly emerging error of closure, 
and of the sudden change of behavior exhibited by the 
ratios fj{t) which is detailed in Fig. [ill These errors 


drive the solution to the Gross-Pitaevskn equation (104) 
off its intended track, to such an extent that any tangible 
connection to the actual Wparticle dynamics appears to 
be lost. 

One might hope that, similar to the effect achieved 
when going from N = 100 in Fig. to = 1000 in 
Fig.0 the loss of t-coherence recorded in Fig. could 
be counteracted by increasing the particle number still 
further. However, this is not the case: Figure il¬ 
lustrates the system’s response when N = 10000, with 


all other parameters equaling those employed in Fig. 10 



FIG. 11: (Golor online) Real and imaginary part of the ratio 
/i for the pulse monitored in Fig. |10[ Observe the change of 
the ordinates’ scale in comparison with Fig. [^ 


FIG. 13: Stiffness |i?i(t)| during the rise of the pulses moni¬ 
tored in Fig. jlOj where N = 1000 (a), and in Fig. 12 where 
N = 10000 (b). Observe the scaling with the particle number. 
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While Fig. testifies that now the stiffness deviates 
from unity by not more than 2 x 10“^ during the rise 
of the pulse, apparently scaling with 1/iV, it then drops 
even more sharply, indicating a “sudden death” of the 
mean field, ie., a dynamically induced sudden loss of 
t-coherence, as corresponding to a sudden destruction 
of the macroscopic wave function [32] ■ With the help 
of the following semiclassical deliberations we will argue 
that this loss of t-coherence persists in the formal limit 
TV —)■ 00 , when a is kept constant: In the situation scru¬ 
tinized in Figs. [To] and [T2| there is no chance to describe 
the system’s time evolution correctly by means of the 
Gross-Pitaevskii equation, not even for arbitrarily large 
particle numbers. 


D. Semiclassical interpretation 

We now exploit the circumstance that the dynamics 
gener ated by the time-dependent Gross-Pitaevskii equa¬ 
tion (1041 are equivalent to those of a driven nonlinear 
pendulum [481149] . Starting from the polar representa¬ 
tion 




Cn = C4 e 


(117) 


of the mean-field amplitudes, and introducing their im¬ 
balance 


^=|C1P-|C2P 

and the relative phase 

(j) = -§2- 'i?l, 


(118) 


(119) 


the equation of motion (1041 readily yields 

■ * * 

1 — = C1C2 - C1C2 


Observing 


where Eq. (122) has been used. On the other hand, one 
deduces 


+C1C2) 

dr 


= 2 


a (|ciP - |c2p)-f/(r) (cic;-CiC2) (125) 


from the Gross-Pitaevskii equation (|104|), having written 

(126) 


fi X M('r) . /W 

/(r) = —s,n(-r 


SO that 

d 


— (ciCt-I- clc 2 ) = —2 az + fir) \/l — sin (j) . (127) 
dr L J 


Gombining Eqs. ( |124| ) and (1271 then leads to 
z 


— = 2 az ■ 
dr 


\/l - 


cos(^ + 2/(r) . 


(128) 


These equations of motion (122) and (128) for the mean- 


field imbalance and the relative phase constitute a pair of 
Hamiltonian equations derived from the classical Hamil¬ 
tonian function 

H-a\p{z, 4>, t) = az^ — v^l — cos (j) + 2 z/(t) ( 129) 

in which z plays the role of a momentum variable, and (j) 
that of its canonically conjugate coordinate, and which 
can thus be interpreted as the Hamiltonian of a nonlin¬ 
ear pendulum with mass propotional to 1/a, and with 
momentum-dependent length, which is driven by the ex¬ 
ternal force /(t): Evidently, one has 


dcj) 

dr 


dH, 


nip 


dz 


dz 

dr 


dH, 


nip 


d(j) 


(130) 


In passing, we remark that for |z| <C 1 this pair reduces 
to 


= -|ci||c2|(e‘^-e-'^) . 

(120) 

dcf) 



dr 



dz 

l-z2=4|cinc2p. 

(121) 

dr 


^ « 2/(r) 


— — sin (p , 


(131) 


this becomes 


dz 

dr 


= — \/l — z^ sin (j) . 


Similarly, one has 

C 1 C 2 + clc 2 = \l\ — z^ cos (j) 


( 122 ) 


(123) 


as corresponding to the equations for the phase and the 
current across a macroscopic superconducting Josephson 
junction m- 

This classical viewpoint underlines the significance of 
the extension of the Lipkin-Meshkov-Glick model (68) 


by the time-dependent drive (72): When /(r) = 0, the 


Hamiltonian (1291 represents a dynamical system with 


giving 

d 


(ciCo -I- cT C 2 ) = z sin (f) cos (j) — \/l — z"^ sin (f) 
dr dr 


d(/ 

~ •) 

( 124 ) 


one single degree of freedom, possessing energy as its 
integral of motion. Adding a time-dependent force is 
tantamount to adding a further degree of freedom not 
accompanied, in general, by a further integral of motion, 
so that the driven nonlinear pendulum is non-integrable 
in the sense of classical Hamiltonian mechanics [551153] , 
giving rise to chaotic motion. This raises the question 
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FIG. 14: (Color online) Grey-scale Husimi snapshots of the evolving A'^-particle state for the pulse monitored in Fig. for 
which N — 100 and /Tmax/fl = 0.9, superimposed onto the instantaneous Poincare surfaces of section (blue) at successive times t 
when — 0.0, 0.6, 0.83, 0.86, 0.88, 0.9, 0.6, 0.0 (upper left to lower right). Here the quantum state initially remains 

associated with an adiabatic invariant corresponding to a closed curve in the regular island, as selected by the Bohr-Sommerfeld 
condition (1351 with n = 0, but then spreads over the stochastic sea when the island becomes so small that the invariant is 
destroyed. 


how the actual, linear 7V-particle system behaves when 
its nonlinear mean-field descendant becomes chaotic. 

A useful link between the A^-particle level and its 
mean-field description in terms of the Hamiltonian sys¬ 
tem ( |130[ ) is provided by the 5'17(2)-coherent states (78): 
Taking the expectation value of the A^-particle Hamilto¬ 
nian (71) with respect to these states, one obtains 




= a \i — 


Nhn/2 

^1 — (cos^ d -I-1) — a/I — COS^ 0 cos (j) 


+2 cos 0 sin(tLit) 


(132) 


which equals the classical Hamiltonian (129) up to an 


irrelevant constant if one neglects terms of order 0{1/N), 
and sets 

z = COS0 . (133) 

Since cosO = cos^(0/2) — sin^(0/2), this is in perfect 
correspondence with the representation (76) of the N- 


fold occupied single-particle state which underlies the 
Af-coheren t sta te \0,(f))N on the one hand, and with the 
definition (1181 on the other. Therefore, the squared pro¬ 
jections 




referred to as Husimi distributions, quantify the degree 
of affinity of the A^-particle state |'I'(t)) with the phase- 
space point (z, V?) at the moment t [iSl |32]. This 
observation leads to the desired semiclassical view on the 
pulse dynamics recorded in Sec. |HI C| We take the so¬ 
lution |'I>(<)) to the Af-particle Schrodinger equation at 
certain moments t, compute the associated Husimi distri¬ 
butions (134), and superimpose these distributions onto 


( 134 ) 


the corresponding classical phase-space portraits, that is, 
onto the Poincare surfaces of section of the accompany¬ 
ing pendulum which is periodically driven with amplitude 
fj,/n “frozen” at the value reached at the moment t under 
study. These surfaces of section are computed by solving 
the classical equations of motion for a representative set 
of initial conditions with fixed driving amplitude /r/H, 
and by plotting the resulting phase-space points strobo- 
scopically after each driving period 2 Trfui. 

Figure [T^ shows a series of snapshots obtained in this 
manner for the pulse previously studied in Fig. for 
which N = 100 and = 0.9. Initially, when 

fi/n = 0, the accompanying classical pendulum is in- 
tegrable, as is reflected by a phase-space portrait which 
is stratified into closed curves remaining invariant under 
the Hamiltonian flow. The initial quantum state consid¬ 
ered here, which is the ground state of the undriven sys¬ 
tem (68) as depicted in Fig.[^ is semiclassically linked by 
means of the WKB-construction to that invariant curve 
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FIG. 15: (Color online) Grey-scale Husimi snapshots of the evolving A'^-particle state for the pulse monitored in Fig. for 
which N = 1000 and /imax/fl = 0.9, superimposed onto the instantaneous Poincare surfaces of section (blue) at successive 
times t when = 0.0, 0.6, 0.83, 0.86, 0.88, 0.9, 0.6, 0.0 (upper left to lower right). While the Poincare sections are the 

same as those displayed in Fig. |14[ now the curve 70 required by the Bohr-Sommerfeld condition ( |135[ ) encircles an area which 
is reduced by a factor of 10 , and thus remains confined in the regular island during the entire pulse, enabling adiabatic following 
of the A'^-particle quantum state. 


7 „ surrounding the central fixed point which is selected 
by the Bohr-Sommerfeld condition m m im [51] 


1 

27r 



— ^eff 



with effective Planck constant 


^eff — 


2 

N 


(135) 


(136) 


and with quantum number n = 0, so that its Husimi 
distribution appears concentrated around that curve 70 . 
When the driving amplitude is increased, a large fraction 
of the invariant curves is destroyed giving way to chaotic 
motion, while those surrounding the central fixed point 
are deformed, but remain preserved if does not be¬ 
come too large, forming a regular island embedded in a 
chaotic sea. If the pulse’s amplitude increases on a time 
scale which is slow compared to the period 27r/a;, these 
preserved curves represent adiabatic invariants to which 
the time-evolving wave function remains tied in a WKB- 
type manner. This is what explains the scenario depicted 
in Fig. |14| For low driving strength the evolving quan¬ 
tum state clinges to its adiabatic invariant still contained 
in the regular island, but with increasing amplitude the 
island becomes so small that the required invariant does 
no longer exist. Then the Husimi projection of the quan¬ 
tum state, having nothing left it can cling to, spills out 


more or less over the entire phase space, leading to a 
final quantum state which contains many eigenstates of 
the undriven junction ( 68 ). 

With this background knowledge the difference be¬ 
tween the pulses recorded in Fig. on the one hand, and 
in Fig.j^on the other, becomes intuitively clear, recalling 
that for the latter pulse the maximum scaled amplitude 
Mmax/^ = 0-9 has been maintained while the particle 
number has been increased to = 1000. Because this 
implies that the effective Planck constant (136) is reduced 
by a factor of 10 , the required phase-space curve 70 now 
encircles a correspondingly smaller area. As shown in 
Fig.[^ it therefore fits into the shrinking and re-growing 
regular island during the entire pulse, albeit just barely 
so at its maximum. Thus, the quantum state evolving 
under the influence of the pulse remains semiclassically 
associated with an adiabatic invariant which is preserved 
from the pulse’s beginning until its end, and thus enables 
adiabatic following. As a result, the final A^-particle state 
here closely resembles the initial one. 


A mere glimpse on Fig.j^then suffices to grasp why no 
further enhancement of the particle number whatsoever, 
together with the implied reduction of the interparticle 
interaction in order to maintain a constant value of a, 
could possibly bring one back to the Gross-Pitaevskii 
regime under the conditions of Figs. [TOj and |T^ Here the 
pulse possesses the maximum amplitude p-max/^^ = 1-5, 
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FIG. 16: (Color online) Grey-scale Husimi snapshots of the evolving A^-particle state for the pulse monitored in Fig. 10 for 
which N = 1000 and /Tmax/fl = 1.5, superimposed onto the instantaneous Poincare surfaces of section (blue) at successive 
times t when = 0.0, 0.73, 1.0, 1.25, 1.43, 1.5, 1.43, 0.0 (upper left to lower right). Here the regular island is completely 

destroyed in the middle of the pulse, implying that there is no adiabatic invariant capable of carrying a macroscopic wave 
function. Thus, even for arbitrarily large particle number the state becomes so complex that a Gross-Pitaevskii description is 
impossible. 


which is so strong that the regular island which still hosts 
the adiabatic invariant at the early stages of the pulse is 
destroyed completely when the envelope approaches its 
maximum. Hence, no matter how large the particle num¬ 
ber N and, consequently, how small the effective Planck 
constant (136), there is no preserved adiabatic invariant 
which could possibly carry a macroscopic wave function. 
Instead, Fig. [T6| visualizes that upon destruction of the 
island the quantum state aquires a degree of complexity 
which strikingly contrasts the simplicity of the initial con¬ 
dition. Obviously this sudden increase of complexity, as 
corresponding to a sudden loss of order, reflects the sud¬ 
den drop of stiffness observed in Figs. [T0| and [T^ signaling 
the sudden death of the mean field. We suggest that this 
scenario is quite general, so that it would be experimen¬ 
tally accessible, for instance, with pulsed Bose-Einstein 
condensates in anharmonic trapping potentials |32j . 


IV. DISCUSSION 

One usually takes recourse to a mean-field theory, 
which in many cases is still tractable by numerical means, 
to obtain information on intricate many-particle dynam¬ 
ics, which in general do not lend themselves to direct sim¬ 
ulations. In our approach to the time-dependent Gross- 
Pitaevskii equation this common strategy has been re¬ 


versed, in order to clarify questions concerning the emer¬ 
gence and possible persistence of an order parameter in 
a system of N Bosons exposed to external forcing: As¬ 
suming that we know its exact many-particle wave func¬ 
tion f['(ri,..., r^v; t), we have asked how this knowledge 
can be used for constructing the associated macroscopic 
wave function, provided that it exists, in the form of 
a mean-field amplitude $(r;t). The answer, formalized 
in Eqs. (33), (38), (39), and (40), is conceptually in¬ 
teresting insofar as it does neither involve the formal 
limit V —> 00 , nor the notion of spontaneous symme¬ 
try breaking: Take the initial V-particle state, annihilate 
a particle, and propagate both the initial state and the 
subsidiary {N — l)-particle states thus obtained in time; 
a candidate mean-field amplitude then is introduced by 
taking the matrix elements of the annihilation operator 
with these evolving states. 


The starting point of this construction resembles the 
specification of a condensate wave function by Lifshitz 
and Pitaevskii [23j, but the additional consideration of 
time evolution in response to external forcing brings 
about a new twist. Namely, in order to decide whether 
or not the candidate actually is a proper mean-field am¬ 
plitude, which is tantamount to the question whether or 
not a macroscopic wave function of the driven system still 
does exist, it is not sufficient to follow solely the trajec¬ 
tory of the given V-particle state in Fock space. Rather, 
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one has to compare this given trajectory to other tra¬ 
jectories initially close to it, differing by one in particle 
number, and to check to what extent these initially close 
trajectories diverge in the course of time. This is what 
is effectively quantified by the projection (65): If the ab¬ 
solute value of this scalar product i?(r; tjequals unity 
to good approximation, the state obtained by first evolv¬ 
ing the initial condition in time and annihilating a par¬ 
ticle at a later moment differs merely by a phase factor 
from the states obtained by annihilating first and evolv¬ 
ing thereafter. This property of t-coherence, indicated by 
|i?(r; t)| = 1, means that initial states differing by one in 
particle number move in some sense parallel to each other 
in Fock space. In this case the flow in Fock space can be 
considered as stiff, similar to laminar flow in fluid me¬ 
chanics. This is the quality of “simplicity” required for 
ensuring the presence of a genuine mean-field amplitude: 
If \R{r; t)| = 1 at least to a good approximation, the can¬ 
didate actually qualifies as a macroscopic wave function 
and obeys the time-dependent Gross-Pitaevskii-equation; 
if not, it has no immediate physical significance. 

The model of the driven Josephson junction studied 
in Sec. m suggests that the occurrence of chaotic solu¬ 
tions to the Gross-Pitaevskii equation reflects the loss of 
t-coherence on the iV-particle level. Our findings thus 
extend previous studies [501155] which have cast doubt 
on the validity of a Gross-Pitaevskii-type mean-field ap¬ 
proach under chaotic conditions. It appears that the 
distinction between order and chaos, which has been 
explored in great depth in classical mechanics jSS] [SO] , 
may have implications of its own in quantum many-body 
physics. In particular, the scenario depicted in Fig. 
indicates in a striking manner that a macroscopic wave 
function can be destroyed almost instantaneously upon 
entering a chaotic regime. 

Recent pioneering works which have taken up the in¬ 


vestigation of many-body quantum chaos have consid¬ 
ered 5-kicked condensates [241 - 129] . In contrast, here we 
have studied the response to a sinusoidal force with a 
“slowly” varying envelope. One of our most important 
results, albeit obtained for one particular model system 
only, consists in the observation that adiabatic follow¬ 
ing to such slowly changing driving forces allows one to 
preserve a pre-existing macroscopic wave function, and 
to transport it almost without loss of t-coherence into 
the regime of strong driving, as shown exemplarily in 
Fig. [7] This finding almost provides a blueprint for gen¬ 
erating Floquet condensates [S2|. More generally, it may 
be of interest for guiding further experiments with time- 
periodically forced Bose-Einstein condensates intended 
to engineer novel systems which may not be accessible 
without such forcing [55H75] : If it is possible to preserve 
maximum stiffness, or t-coherence, even in the presence 
of strong forcing, it should also be possible to actively 
manipulate macroscopic wave functions by applying suit¬ 
able coherent control techniques. The evidence collected 
in the present work clearly indicates that this road is 
viable. 
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